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Abstract. The coexistence of multiple types of interactions within social, 
technological and biological networks has moved the focus of the physics of complex 
systems towards a multiplex description of the interactions between their constituents. 
This novel approach has unveiled that the multiplex nature of complex systems has 
strong influence in the emergence of collective states and their critical properties. Here 
we address an important issue that is intrinsic to the coexistence of multiple means of 
interactions within a network: their competition. To this aim, we study a two-layer 
multiplex in which the activity of users can be localized in each of the layer or shared 
between them, favoring that neighboring nodes within a layer focus their activity on 
the same layer. This framework mimics the coexistence and competition of multiple 
communication channels, in a way that the prevalence of a particular communication 
platform emerges as a result of the localization of users activity in one single interaction 
layer. Our results indicate that there is a transition from localization (use of a preferred 
layer) to delocalization (combined usage of both layers) and that the prevalence of a 
particular layer (in the localized state) depends on their structural properties. 
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1. Introduction 

During the last 15 years many statistical physics methods and nonlinear models have 
suffered a reformulation in order to take into account non-regular interaction patterns 
[DEI [3]. This reformulation is rooted in the availability of datasets capturing the 
relationships among the constituents of macroscopic systems of diverse nature, such 
as technological, biological and social ones, and their successful description in terms of 
complex networks As a result, many important collective phenomena taking place 

in these systems, such as synchronization [6], epidemics [7], cooperation and consensus 
[8] among others, have been revisited under the paradigm of complex networks [HI fin] . 

As data gathering techniques increase their resolution new properties of the 
interaction patterns in complex systems are captured. Main features include the spatial, 
temporal and multiplex nature of interaction networks. This latter ingredient has greatly 
focused the attention of network science in the last years leading to a number of works 
about the structure and dynamics of multilayer and multiplex networks mm- 

Multiplex networks [I3] are often described as the framework for capturing the 
interactions among a set of elements (nodes) when these interactions can take different 
forms or be established through different means, each of them defining a network often 
referred to as interaction layer. Thus, a multiplex network can be seen as a collection 
of layers so that each node is represented in all of them (see. Fig. d]). Typical examples 
of multiplex networks are transportation systems [H [ig [m [m, in which different 
transportation modes can be used to travel between cities or urban areas, or social 
systems [18] , in which individuals can choose among different means and communication 
platforms for interacting with each other. 

Most of the times the different interaction layers forming the multiplex coexist 
in a competitive way [m [201 [H]. For instance think of two of the most important 
applications for mobile communication, such as WhatsApp and Line, here the 
competition relies on the usage of each platforms. The more users decide to use one 
platform the more value has the platform. From the point of view of users the final 
choice between these two platforms relies on two main issues. Obviously, the intrinsic 
quality of the platform plays a key role in the final decision of individuals. However, 
there is social added value that comes from the degree of usage of the platform among 
the acquaintances of an individual. It is thus interesting how this local and context- 
driven decisions affect the onset of a collective state, here represented as the localization 
of the multiplex activity in one of its layers. 

The article is organized as follows. In Sec. |2] we present the mathematical 
formulation of the multiplex network model and derive the relevant equations to be 
solved. In Sec. [3] we show the numerical results corresponding to the competition 
between different interaction layers. Finally, in Sec. ID we round off the manuscript 
with the conclusions of our work and the future perspectives that it may open. 
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Figure 1. Schematic representation of a two-layers multiplex network. The multiplex 
is composed of 6 nodes so that each node appears in each layer and it is connected 
(dashed lines) with its representation in the other layer. When layers compete, each 
user must choose the layer in which is active (black), thus remaining inactive (white) 
in the other one. 

2. Mathematical formulation of multiplex layers competition 

The physical framework used to study the layer-layer competition relies on a two-states 
model (reminiscent of an Ising-like model but much simpler). In this case, nodes are 
two-states (up and down) systems and the interaction neighborhood of a node depends 
on its (up or down) state. Thus, at the macroscopic scale, this translates into the 
existence of two interaction layers, one associated with the up state and another with 
the down one. In the following we hrst describe in 12. li the general Hamiltonian capturing 
the interactions between the nodes of a multiplex of L layers to particularize in 12.21 to 
the case of competitive interactions in a two-layers multiplex. 

2.1. General formulation: Multiplexes of L layers 

In general, a multiplex network is composed of L network layers of N nodes each. Since 
each individual i is represented in each of the L networks, each pair of networks a and f3 
are interconnected by N links connecting the pair of nodes that represent the same 
individual i. This setup can be seen as a collection of independent layer or platforms 
available for the communication between individuals (such as WhatsApp, Line, Tango, 
etc) or for data-sharing (such as Dropbox, iCloud, Box, etc), being the nodes of the 
layers, the users, and the links within a layer the connections established by the users 
via a particular platform. 

The availability of different platforms oriented towards the same goal poses a natural 
competition for the choice of the users. The essential ingredients of this competition 
can be casted in a mathematical formulation in which the state of a node i in layer a 
can be explained as the probability that node i is active in layer a (for communication 
platforms) or the fraction of resources that node i share with its neighbors in this 
layer (for data-allocation systems). In this way, the state vector of a node i in the 
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interconnected multilayer network is denoted by pi = {p[^\ ■ ■ ■ ,Pi^'^), together with the 
constraints that the sum of the probabilities of hnding a node active in each layer is 
equal to 1: 

L 

= z = (1) 

The state of the multiplex can be represented by a L x iV matrix P = 
[(pi)!,..., {pnV] ■ The state matrix, accounting for the intensity of interaction between 
the nodes is given by all possible products between the states, he., by the LN x LN 
matrix S = P $§ with elements = p^f^p^P ■ In addition, we define J as a 

LN X LN interaction matrix capturing both intra-layer and inter-layer links: 

J = + D(8)I , (2) 

where I indicates the N x N identity matrix, D is the L x L matrix accounting for the 
network of layers ra and Wp is the weight of the interaction between nodes i and j 
in layer a. More specihcally, this topology describes nodes that are present in multiple 
layers simultaneously and inter-layer connections are allowed only between a node and 
its counterparts in the other layers. If a and are indices indicating two given layers, 
the block matrix structure of J can be indexed by four indices, two for nodes and two 
for layers, he., by ■ 

Therefore, given the ensemble of all possible states, he., the set {P} of all matrices 
satisfying constraints Eq. o, we can define the Hamiltonian of a specihc conhguration 
P as: 

L N 

w(p)=- E E . (3) 

Q;,/3=1 ij = l 


2.2. Competition in 2-layers multiplexes 

For the sake of simplicity, hereafter we consider the case of the competition in a multiplex 
composed of two layers, i.e., L = 2. In this specihc case, following Eq. ([I]), the state 
of a node is completely determined by its probability of being active in one of the two 
layers, e.g., the hrst one. Thus, the state of the whole multiplex can be described by the 
vector p = (pi,... ,pn)-i where for simplicity we have omitted the layer index explicitly. 
Moreover, we also consider a uniform and undirected connection between the two layers, 
so that the interaction strength of a node i in layer 1 and its counterpart in layer 2 is 
captured by the parameter J^. 

To incorporate the competition between layers we consider two essential ingredients 
of the interactions at the local level. On one hand, the communication between two 
agents that are connected within one of the layers is more efficient when both of them 
are always active in this layer or when they allocate all of their shared resources in the 
same platform. However, since the sets of contacts an individual has in the two layers 
are, in principle, different, by splitting the activity between the two layers an individual 
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will increase the number of simultaneous contacts. The Hamiltonian capturing these 
two ingredients can be obtained from the general one in Eq. (|3]) as: 

N N N 

H{p) = - Wi - E W,f{l-p.){l-n) - -k).(4) 

i,j=l *j'=l *=1 

From this Hamiltonian it becomes clear that the first two terms in the right are 
those favoring the localization of the activity of each individual in layer one and two 
respectively. In its turn, the third term favors the splitting of the activity of each 
individual. 

The relative importance of this third term with respect to those favoring the 
localization of the activity within a single layer is controlled by the inter-layer coupling 
Jx- Note that the limit 3> 1 means that nodes are prone to combine their activity 
in both layers which, for instance, in the case of data sharing or mobile communication 
platforms would represent information (pictures, files, tweets, etc) that can be easily 
transferred from one platform to the other one. On the contrary, the case 1 

implies that a simultaneous use of platforms is hard to achieve. 

3. Results 

Having introduced the mathematical framework, our goal is to study the competition 
between the two layers as a function of the inter-layer strength and the structural 
patterns of each of the network layers. To this aim, it is useful to check the behavior in 
the two asymptotic limits: 3> 1 and = 0. First, when S> 1 the first two terms in 

the Hamiltonian Eq. (jlj) become negligible, so that the configuration of minimum energy 
is achieved for pi = 1/2 Vi, i.e., when the individuals split their activity between the two 
layers. On the other hand, for Jj, = 0 the multiplex becomes a set of two independent 
networks and the configurations localized in the first layer {p = 1, i. e. pi = 1 Vi) and the 
second one (p = 0) compete. In this case the minimum energy conhguration is achieved 
by concentrating all the activity in the layer a with the largest total strength: 

N 

,(«) = Y, . ( 5 ) 

i,i=l 

For the particular case of unweighted networks this means that the layer with the largest 
average degree (or largest number of links) will concentrate all the activity when = 0. 
From now on, we will consider that so that in the absence of inter-layer 

interactions the activity focuses on the first layer: p = 1. 

Considering these two asymptotic behaviors we are thus interested in characterizing 
the transition from the localized activity regime at small values of inter-layer coupling 
Jx to that of mixed one for large Jx- A first proxy is to check when the fully localized 
solution {e.g. p = 1) ceases to be the one with the minimum energy H = —(a = 1 in 
the case of p = 1). To this aim, we calculate the gradient of the multivariate Hamiltonian 
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H{p): 


dH 

dpi 


N N 

-2 K’Pi + 2 E ’(1 - ft) - 2i(l - 2 b) . 

i=i i=i 


( 6 ) 


The values of these derivatives for the localized solution in the hrst layer {p = 1) become: 


dH 


dpi 


p=i 



(7) 


where ^ is the strength of node i in layer a. For values of smaller than 

i=i 

the strength of all nodes in the hrst layer, the derivatives at p = 1 are all negative 
and thus the gradient points to the interior of the hypercube [0,1]^, which contains all 
the possible feasible states of the system. This means that the energy of the system 
around p = 1 is always increased for any small change of p inside the hypercube, showing 
that p = 1 has minimum energy whenever the inter-layer coupling is below its critical 
value: 


U rr. 


min 


mm * 


( 8 ) 


For unweighted networks it reduces to the minimum degree of the nodes in the hrst 
layer, Above this critical inter-layer coupling the minimum energy moves from 

p = 1 to a new position inside the hypercube [0,1]^, thus starting to distribute the 
activity between the two layers. 

From a mathematical point of view, the hnding of the state with minimal energy is a 
quadratic (the Hamiltonian) programming problem with linear equality (normalization 
of the probabilities) and inequality (probabilities in range [0,1]) constraints. In general, 
for two layers, the candidate minimum of the Hamiltonian is calculated by setting 
dH/dpi = 0 Vi, which can be expressed as the following linear system: 

[2JJ-(wW + W(2))]p= 41-^2) _ ^9^ 


However, its solution fL does not always fulhll the constraints, does not constitute a 
minimum, or even both conditions fail at the same time. When any of these happens, 
the minimum is placed in the boundaries of the [0,1]^ hypercube, with at least one 
probability equal to 1 or 0. Supposing p* is inside the hypercube, it is a minimum if the 
Hessian matrix, with components 

II-= 2 (2J...,-<>-<)) , (10) 

is positive dehnite. When the Hessian is not positive dehnite, the quadratic 
programming problem becomes NP-hard, and no polynomial time algorithm is known 
to solve it. Since the Hessian is proportional to the matrix of the system in Eq. ([9]), a 
positive dehnite Hessian implies the system has a non-singular matrix and thus a unique 
solution. 

Summarizing, we proceed as follows to hnd the state with minimum energy. For 
each value of the coupling Jx, we hrst solve Eq. (I9]) and obtain a solution jp. Then, we 
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Figure 2. Magnetization for multiplex networks consisting of two Erdos-Renyi (ER) 
layers and varying mean degree. All ER networks have 200 nodes and minimum degree 
kmin = 1- The red circles indicate the minimum values of Jx so that minimum energy 
configurations can be calculated from Eq. For values of Jx to the left of the red 
circle, PSO has been used to compute the solutions. In the inset, zoom to show when 
activity starts to be distributed in both layers. 


check ii fL G [0,1]^ and if all the eigenvalues of the Hessian Eq. flTOl) are positive. If 
both conditions are fulhlled, pL is the ground state for this value of the coupling and we 
have hnished. Otherwise, a heuristics is needed to solve the problem. We have chosen 
Particle Swarm Optimization [23] for its simplicity, ability to cope with continuous and 
bounded variables, and outstanding performance in many helds |24] . 

In order to represent the ground state of the Hamiltonian for each value of Jx we 
make use of magnetization M{p): 

1 ^ 

( 11 ) 

i=l 

that characterizes the level of activity between layers: M = 1 when all the activity is 
concentrated in the hrst layer [p = 1) and M = — 1 when it is concentrated in the second 
layer (p = 0). Obviously, when activity is shared between the two layers {p = 0.5) the 
magnetization vanishes, M = 0. 

In Fig. [2]we show the magnetization as a function of Jx for three multiplex networks 
composed of two Erdos-Renyi (ER) layers. The layers of the three multiplex networks 
have average degree {k) = 4, 6 and 8 respectively. The layers were produced by means 
of the algorithm introduced in [25] that allows to interpolate between ER and scale-free 
networks and, as a byproduct of the procedure, enables to control the minimum degree 
of the resulting graphs. In this way we have set the minimum degree to /cmin = 1 in the 
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Figure 3. Magnetization for multiplex networks consisting of two Erdos-Renyi (ER) 
layers and varying mean and minimum degrees. All ER networks have 200 nodes. The 
red circles indicate the minimum values of Jx so that minimum energy configurations 
can be calculated from Eq. (I9|). For values of Jx to the left of the red circle, PSO has 
been used to compute the solutions. In the inset, zoom to show when activity starts 
to be distributed in both layers. 


three cases. However, the method in [25] produces (for a given value of kram) networks 
of identical strength. Thus, we take the hrst network layer an add 0.05 x N links at 
random to ensure that The three curves in Fig. [2] display the transition from 

localized activity for small values to mixed one for large ones. Interestingly, the inset 
shows that localized activity is lost as soon as = 1 = fc min in agreement with our 
former estimation. From this point the decay of M is slower for those multiplexes with 
larger mean degree {k) in the layers. 

These results are further corroborated in Fig. [3l In this case we show again three 
two-layers multiplexes (again composed of coupled ER-ER networks) with the same 
average degrees as in Fig. [2] {{k) = 4, 6 and 8) but with different values of the minimum 
degree, fc min = 1, 2 and 4 respectively. This latter feature is revealed in the inset of the 
plot where we show that the state of localized activity (M = 1) is no longer the ground 
state of the multiplex as soon as > /Cmin = 2 and 4. On the other hand, as in 

Fig. O the hnal mixed activity state is achieved hrst for those multiplexes with smaller 
average degree. 

In both Figs. [2] and [3] the insets show that the evolution of M{Jx) close to M = 1 
shows different cusps. In particular, these cusps appear at integer values of Jx, which at 
the same time correspond to the values at which the number of nodes for which Eq. fl71) 
changes sign increases. The effect is similar to that at kmm- many nodes have pi = 1 
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until Jx is equal to ki in the first layer. This rule is exact only for nodes with ki = femin, 
but also holds for an important fraction of nodes when ki > /cmin- The collective effect 
is reflected in the form of cusps in the curve M{Jx). 

4. Conclusion 

In this work we have introduced a model to analyze how layers compete for the activity 
of users in a multiplex network inspired in the simultaneous interplay of communication 
and data-sharing online platforms. To this aim, we have focused on a multiplex 
composed of two layers and we have relied on a two-states model in which each of 
the two states of a node are associated to be active in the top layer and the bottom 
one respectively. At variance with the usual Ising model in a network |26], here a node 
interacts only with those neighbors in the layer it is active. We have set two competing 
mechanisms, one favoring activity localization and another favoring the splitting of the 
node’s activity between the two layers. On one hand, if a node and all its neighbors 
in a layer are active in it, this would favor an efficient communication between them. 
However, when a node splits its activity between the two layers this would favor the 
passage of information from a neighbor in one layer to a neighbor in the other one, 
thus maximizing the outreach of information. The competition between these two 
mechanisms is controlled by the inter-layer coupling which can be seen as the ability 
that node has to pass information from one layer to the the other. 

Our results show that, regardless of the average connectivity and total strength of 
the leading layer (the one that focuses all the activity for small inter-layer coupling) it 
is its minimum degree what causes that nodes start to use the other layer, i.e. controls 
the onset of the transition from localized to mixed activity. On the other hand, it is 
the average degree of the leading layer what controls when the state of full mixing is 
reached. These two results point out that the transition from localized to mixed activity 
occurs via a cascade from poorly connected nodes in the leading layers (the ones that 
obtain more benefits from leaving first the leading layer) to those highly connected ones 
(being the ones that are less prone to leave the layer in which they are well-connected). 
Thus, the larger the average degree of the leading layer, the more inter-layer coupling 
is needed to persuade all the degree classes to leave the localized state. 

We expect that the simple model introduced here will stimulate more research about 
the coexistence and competition of interaction layers in multiplex networks making 
possible the characterization of how and when the coexistence of different layers in real 
multiplex systems is possible. Future research avenues include the study of other types 
of layer topologies and the presence of correlations between the degrees of a node in 
different layers. Moreover, a more challenging problem is the competition in multiplexes 
composed of more than two layers, characterized by the Hamiltonian in Eq. ([3]). It is 
clear, that the existence of multiple parameters for the interaction between the L layers 
poses a mathematical and computational difficulty. On the other hand, this general 
framework provides with the interesting scenario in which many different transitions 
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between the localization in different layers are observed due to the multiple competition 
between them. It is also interesting to note that we have tackled the analysis of the two- 
state model by considering continuous variables, {pi}, associated to each node. However, 
another possibility is to consider binary states for the nodes, as in Ising-like models, so 
that metastable states can show up due to the multi-stable character of the node states 
[271128]. 
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